Quantum trajectory approach to circuit QED: Quantum jumps and the Zeno effect 
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We present a theoretical study of a superconducting charge qubit dispersively coupled to a trans- 
mission line resonator. Starting from a master equation description of this coupled system and using 
a polaron transformation, we obtain an exact effective master equation for the qubit. We then use 
quantum trajectory theory to investigate the measurement of the qubit by continuous homodyne 
measurement of the resonator out-field. Using the same porlaron transformation, a stochastic mas- 
ter equation for the conditional state of the qubit is obtained. From this result, various definitions 
of the measurement time are studied. Furthermore, we find that in the limit of strong homodyne 
measurement, typical quantum trajectories for the qubit exhibit a crossover from diffusive to jump- 
like behavior. Finally, in the presence of Rabi drive on the qubit, the qubit dynamics is shown to 
exhibit quantum Zeno behavior. 

PACS numbers: 03.65.Yz, 42.50.Pq, 42.50.Lc, 74.50.-|-r, 03.65.Ta 



I. INTRODUCTION 

Continuous- in-time measurement theory [1, 2], or 
quantum trajectory theory, describes how an observer's 
state of knowledge of a quantum system (known as the 
conditional state) evolves given a measurement record. 
Even in the absence of classical noise, such a trajectory 
follows a stochastic path in time, with the randomness 
being due to quantum uncertainty. These stochastic tra- 
jectories are either diffusive or jump-like in nature. Dif- 
fusive trajectories usually arise when the observable be- 
ing measured is only weakly coupled to the detector [3], 
whereas jump-like behavior occurs when there is a large 
sudden change in the observers knowledge of the system 
state, a typical example of the latter being detection of 
a photon with a photomultiplier [1, 4]. The evolution 
equation for this trajectory is called a stochastic master 
equation (SME) [1, 2, 3, 4, 5]. 

In this paper, we consider measurement in circuit quan- 
tum electrodynamics (QED) [6, 7, 8, 9, 10]. This system 
consists of a cooper pair box, playing the role of an ar- 
tificial atom, dispersively coupled to a 1-D transmission 
line resonator and is the circuit equivalent of cavity QED. 
It has the advantage that the qubit can be fixed at an 
antinode of the resonator which, due to its ID configura- 
tion, has a very large vacuum electric field. This leads 
to very strong Jaynes-Cummings type coupling between 
the qubit and the resonator. This allows to probe a new 
regime of parameter space, where the resonator and qubit 
are strongly coupled via a dispersive interaction and no 
energy is exchanged between them. In this regime, the 
qubit causes a large state-dependent shift of the resonator 
frequency and thus by monitoring the signal transmit- 
ted through the resonator we can infer the qubit's state. 
An important feature of the quasi-one-dimensional circuit 
resonator is that it permits dispersive couplings ^ 10^ 



times larger than for ordinary three-dimensional cavi- 
ties. Coupling of superconducting charge qubits to 3D 
microwave cavities has also been investigated theoreti- 
cally [11]. 

To go beyond previous theoretical works on this sys- 
tem, we employ quantum trajectory theory [1, 2, 3, 4, 5] 
and a polaron transformation [12] to derive an effective 
SME for the qubit. It is well know in quantum optics that 
if the voltage of a cavity is weakly monitored, then the 
evolution of the conditional state of the cavity is given by 
the homodyne SME [1, 13]. Using this SME as a start- 
ing point, we show that if the rate at which information 
is coming out of the resonator is much larger then the 
rate at which information is being lost into unmonitored 
baths, then we can use a polaron transformation to elim- 
inate the resonator and obtain a SME for the qubit only. 

The effective SME which is derived corresponds to a 
weak dispersive measurement of the qubit observable ■ 
This is akin to SMEs derived in Refs. [14, 15, 16, 17, 18, 
19, 20, 21] for a quantum dot monitored by a quantum 
point contact. In Refs. [14, 15, 16, 19], the SMEs are re- 
ferred to as the quantum Bayesian equations. They can 
be seen as the Stratonovich version of the the Ito SMEs 
presented here and in Refs. [18, 21]. We note that one 
can derive a similar equation by adiabatically eliminating 
the resonator from the qubit-resonator master equation 
using a similar technique as presented in [22]. However, 
by using the present method, we are able to derive correc- 
tions to the various system rates. These rates are shown 
to agree very well with the solution of the total resonator- 
qubit conditional state found by numerically solving the 
homodyne SME. 

The paper is organized as follows. The next section 
contains a brief discussion of the circuit QED hamilto- 
nian and of the dispersive approximation. In Sec. Ill, we 
use a polaron transformation to eliminate the resonator 
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from the resonator-qubit master equation and in this way 
obtain an exact master equation for the qubit only. This 
exact effective master equation shows how the resonator 
induces an additional dephasing channel on the qubit, 
whose strength depends on the amplitude of input drives 
on the resonator. This dephasing rate corresponds to 
measurement-induced dephasing found in Ref. [10]. It 
contains the number splitting predicted in Ref. [10, 38] 
and experimentally reported in Ref. [24]. In Sec. IV we 
apply the same polaron transformation to derive an ef- 
fective SME for the qubit. This equation corresponds to 
a weak measurement of the qubit operator with rate 
Fci and extra non-Heisenberg backaction which causes 
random rotations of the qubit around the ct^ with rate 
Fba- Using these results, in Sec. V we investigate the sit- 
uation where the measurement is quantum limited and 
how a measurement time can be defined. In Sec. VI, we 
describe the emergence of quantum jumps in the condi- 
tional state as the measurement strength is increased. In 
particular, we investigate how the spontaneous relaxation 
into the qubit bath reveals itself as a jump when the mea- 
surement is strong, rather than a diffusive decay when the 
measurement is weak. Finally, in Sec. VII, we include a 
qubit control drive in our description of the system and 
investigate how the Zcno effect [25, 26, 27, 28] can be 
observed in this system. We summarize our conclusions 
in Sec. VIII. 



II. CAVITY QED WITH SUPERCONDUCTING 
CIRCUITS 

We consider a Cooper pair box capacitively coupled 
to a transmission line resonator acting as a simple har- 
monic oscillator. This system, illustrated schematically 
in Fig. 1, was first introduced in Ref. [6] and experimen- 
tally studied in Refs. [7, 8, 9, 24, 29, 30]. Measurement- 
induced dephasing was theoretically studied in detail in 
Ref. [10] and the applications to quantum information 
processing investigated in Ref. [31]. 

As described in the above references, the system's 
Hamiltonian in the presence of a microwave drive of am- 
plitude ed{t) and frequency ujd can be written as [6] 

H = — -a, + Tujjro) a + Tiq (a) + a(7+) 

2 ^ ^' (2.1) 

+ h [erf(t)a^e-'"''* + e*ait)ae'^^'] . 

In this expression, ujr is the resonator frequency, iUa the 
qubit transition frequency and g the resonator-qubit cou- 
pling strength. Depending on its frequency, the drive can 
correspond either to measurement of the qubit or can be 
used to coherently control its state. In the dispersive 
regime, when |A| = |wa — Wrj 3> \g\, the effective Hamil- 
tonian Eq. (2.1) can be approximated by [6] 

HcS ~ —^az + hAra^a+hxa' aaz + h [edit)a^ -f e^(t)a] , 

(2.2) 




FIG. 1: (Color online) Schematic layout and lumped element 
version of the circuit QED implementation. A superconduct- 
ing charge qubit (green) is fabricated inside a superconducting 
ID transmission line resonator (blue). Here, we have taken 
the input (LHS) and output (RHS) capacitance of the res- 
onator to be unequal. With the output capacitance much 
larger than the input one, a larger fraction of the photons in 
the resonator will escape on the output side. In this way, the 
signal to be measured is enhanced. 

where — ujr ~ ujd and wc have assumed that the drive 
is far off resonance from the qubit transition frequency 
and moved the resonator to a frame rotating at frequency 
LUd- In this expression, we have defined x = 9^ / ^ as the 
dispersive coupling strength between the resonator pho- 
ton number and the qubit. Moreover, we have taken 
Coa = + X as the Lamb shifted qubit transition fre- 
quency and defined A^. = oj^ ^ oJd- 

It is important to note that the dispersive approxi- 
mation breaks down as the number of photons in the 
resonator approaches the critical photon number Wcrit = 
A^/4g^ [6]. In the present paper, we will work at mod- 
erately low photon number and assume the dispersive 
approximation to hold. That is, we will assume that 
Eq. (2.2) is a valid description of the system. The results 
obtained will therefore also apply to the circuit QED im- 
plementation that uses the so-called transmon qubit [32] . 
A future publication will explore the break down of the 
dispersive approximation [33] at small detuning and large 
photon number. 

III. MASTER EQUATION 

A. Master equation for the combined system 

In the Born-Markov approximation, the master equa- 
tion describing circuit QED takes the usual Lindblad 
form [1, 34] 

m = - ^ [i?eff, Q{t)] + KV[a]Q{t) + ^iV[a]Q{t) 
=CtotQ{t), 

where Q{t) is the state matrix for both the qubit and the 
resonator and 'D[A\ is the damping superoperator defined 
by the mapping 

V[A]q^ AqA^ - qA^A/2. (3.2) 
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The three damping channels are photon loss through the 
resonator (k), qubit decay (71) and dephasing of the 
qubit (7^) [52]. 

There are two distinct contributions to the resonator 
damping: loss of photon from the input port described by 
the rate Kjn and loss of photons at the output port Kout- 
The sum of these two rates k — Kjn + Kout is what appears 
in Eq. (3.1). The advantage of distinguishing these two 
contributions is that, in current experiments, only the 
photons leaking out of the output port are monitored. 
Hence, as we will see later, Kout is related to the rate at 
which we acquire information about the qubit state. 



B. Master equation for the qubit 

In Ref. [10] the solution of the above master equation, 
neglecting energy loss due to 71 (but keeping its effect on 
dephasing), was found to be 

Q{i)= E c„-(t)|i)(j|®|a^(t))(aj(t)|, (3.3) 

where the indices g and e label the qubit ground and 
excited states, respectively. In this expression, the coeffi- 
cients Ci J (t) are given by 

Ce,e(i) = Ce,e(0) 
Cs,g(0 = Cg,3(0) 

Ce g(0)e^''('^"~*''2^*e~*^'<^-''i* "=(s)"s (3.4) 



{ag{t)\ae{t)) 



and |ae(g)(t)) are coherent states of the resonator with 
amplitudes determined by 



ae{t) = -ied{t) - i{Ar + x)o:e{t) - Kae{t)/2 
o^g{t) = -i£d{t) - i{Ar - x)ag{t) - Kag{t)/2. 



(3.5) 



In the expression for Ce^g{t), we see that in addition to 
decay due to dephasing 72 =71/2 + 70, the off-diagonal 
element of the qubit density matrix decays at a rate that 
depends on the field amplitudes Ug and ag. Since the co- 
herent states with these amplitudes act as pointer states 
in the measurement of the qubit, this decay can be in- 
terpreted as measurement-induced dephasing and will de- 
pend on how distinguishable the states \ae(t)) and \ag{t)) 
are [10]. Following Ref. [10], the effect of qubit relaxation 
in this model has been described by Bonzom et al. [35]. 

How these coherent states act as pointer states for the 
qubit is illustrated in Fig. 2, where a phase-space repre- 
sentation of the resonator state Trqubit[£'(i)] is plotted for 
various qubit-field interaction times. In these plots, the 
initial state is taken to be the vacuum for the resonator 
and the symmetric superposition (|e) -I- \g))/\/2 for the 
qubit. A continuous coherent drive is turned on at time 
t = 50 ns to build up photon population of the resonator. 
The three panels of Fig. 2 illustrate how the two coherent 
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FIG. 2: (Color online) In-phase / = Re[{a)] = {a + a^)/2 and 
quadrature Q — Im[{a)] — {ia^ — ia) /2 component of the Q- 
function solution for the resonator state Trqubit[ps]. Here the 
parameters are K/27r = 10 MHz, x/27r = 5 MHz, e/27r = 20 
MHz, and Ar = 0. The measurement drive is turned on at 
t —50 ns and the initial state is jO) (g) (je) + \g))/\/2. For 
illustration purposes, 71 was taken to be zero. A) t = 0, B) 
t = 75 ns, C) Steady-state solution. 



states \ae{t)) and \ag(t)) eventually become separated in 
phase-space. Homodyne detection of the resonator field, 
with the proper choice of local oscillator phase (j), can 
then be used to distinguish between these two coherent 
states and thus readout the state of the qubit. In the last 
panel of Fig. 2 we have introduced the distance between 



the states \ae{t)) and \ag{t)) 
I3{t)=ae{t) 



a, 



Another useful quantity is the angle 
00 = arg(/3). 



(3.6) 



(3.7) 



For example, for the case of Fig. 2C), this angle is Op = tt, 
which corresponds to the quadrature containing the most 
information about the state of the qubit. Indeed, for 
this particular case, all of the information is in the in- 
phase component /, with no information stored in the 
quadrature Q component. 
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Our goal in the remainder of this section is to obtain 
an effective master equation for the qubit by ehminat- 
ing the resonator degree of freedom from the full master 
equation, Eq. (3.1). To achieve this, we first go to a frame 
defined by the transformation 

P{t) = neZ?K(t)] + IlgD[ag{t)] (3.8) 

with D[a\ the displacement operator of the resonator, 

D[a] ~ cxp[Q;a^ — a*a] (3.9) 

and Hj = \j){j\ projectors on the ground and excited 
states of the qubit. This is similar to the polaron trans- 
formation which has been used extensively in various sys- 
tems [12, 36]. For example, this was used in Ref. [37] to 
study a charge qubit coupled to a mechanical oscillator 
beyond the rotating wave approximation. However, here 
we use the transformation Eq. (3.8) on the dispersive 
Hamiltonian, not on the full Jayncs-Cummings Hamilto- 
nian. 

As is shown in Appendix A, applying this transforma- 
tion on the master equation Eq. (3.1) and tracing over 
the resonator state yields the laboratory frame reduced 
qubit master equation 



m 



+ [70 
--Cp{t), 



(3.10) 



where p(t) = TT,-cs[g{t)]. This expression is the main 
result of this section. It is important to note that, within 
the dispersive approximation, this result is exact. 

In this reduced description, we see that the effect of 
coupling to the resonator translates into an additional 
dephasing rate Td{t) given by 



(3.11) 



In addition to dephasing, the qubit transition frequency 
is also modified by the photon population of the res- 
onator. The shifted qubit frequency is given by 



■ UJa 



Bit) 



with 



B{t) ^ 2xRe[ag{t)a:{t)]. 



(3.12) 



(3.13) 



This last term gives rise to the ac-stark shift experi- 
mentally measured in Ref. [9]. In the situation where 
X 2> K, the above leads to the number splitting predicted 
in Ref. [10, 38] and experimentally observed in Ref. [24]. 

Although the above results are analytically exact, we 
have compared the numerical solution of the qubit 's dy- 
namics obtained from the full master equation Eq. (3.1) 
to that obtained from the reduced model Eq. (3.10). This 
is shown in Fig. 3 where the elements of the Bloch vector 
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FIG. 3: (Color online) A) Numerical solution of the three 
components of the Bloch vector as obtained from the full 
master equation: x red dashed, y green dashed-dotted and z 
solid blue. For this simulation, a measurement drive tuned 
at the resonator frequency (A^ — 0) and with envelope 
£d{t) = v4tanh[(t — ton)/cr] is turned on at ton = 20 ns (green 
dashed-dotted hne in panel B). The turn-on time is indicated 
with the vertical dashed hne. B) Purity p (solid blue line) 
and the trace distance d-rr (red dashed line) with respect to 
the effective model derived here. The latter shows the excel- 
lent agreement between the model and full numerical results. 
For these results, we have taken A/2tt = 5 MHz, a — 5 ns, 
Ti — 7 i-is and T2 — 500 ns. The other parameters are the 
same as in Fig. 2. 



are plotted as a function of time. It is useful to note that, 
in terms of the Bloch vector 

p{t) = ^ (1 + x{t)a, + y{t)ay + z{t)cj,) , (3.14) 

where x{t) = Tv[axp{i)\, y{t) = Tr[crj,/?(t)], and z{t) ~ 
Tr[cr2p(t)], the effective master equation, Eq. (3.10), re- 
duces to the simple Bloch equations 



X{t) = -Wac(Oy(i) - [72 + TA{t)]x{t) 

y{t)=^UtXt)-b2 + TA{t)]y{t) 
i(0 --7i(^(0 + l)- 



(3.15) 



Not surprisingly, as seen from Fig. 3, the agreement be- 
tween numerical integration of these Bloch equations and 
full numerical integration of the full master equation is 
excellent. The agreement between these two results can 
be quantified by calculating the trace distance dTr(i) be- 
tween the full and the effective model. The trace distance 
in this case is defined as [39] 

dT^.(t)-iTr[|p(i)-Tr,es[e(0]|], (3-16) 

and is found to be zero at all integration times, up to nu- 
merical round off (truncation) errors. This was checked 
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for a wide range of measurement amplitudes Ed and the 
trace distance was also found to be zero for all verified 
amplitudes. Also in this figure is shown the purity (full 
blue line) 

p{t) = Tv[p[tf] = i[l + x{tf + y[tf + z[tn (3.17) 

In this master equation description, the purity tends 
to 1/2, corresponding to a completely mixed state, due 
to dephasing. On top of this dynamics, relaxation is 
taking the purity to 1 (a pure state) at a rate 71 . For the 
parameters chosen in Fig. 3 relaxation is much weaker 
than (measurement-induced) dephasing, such that only 
the effect of the latter is seen. 

To summarise this section, we have obtained an exact 
master equation description of the qubit dynamics that 
only involves classical solution of the resonator field. In 
addition to giving direct insights in the qubits dynam- 
ics (e.g. measurement induced dephasing and ac-Stark 
shift) this also provides a tool to significantly reduce the 
complexity of numerical calculations. Indeed, only a 2- 
dimensional Hilbert space is now required. We note that 
this effective model was used to analyze and reproduce, 
with exceptional agreement, the experimental results re- 
ported in Ref. [29]. 



IV. STOCHASTIC MASTER EQUATION 

A. Stochastic master equation for the combined 
system 

In this section, we review homodyne measurement of 
the field emitted from the resonator. Using these results, 
we will in the next subsection use the transformation of 
Eq. (3.8) to obtain an effective SME for the qubit only. 

In a given quantum system, if all the decoherence and 
decay channels can be monitored continuously, it is pos- 
sible to describe the system conditioned on the results of 
this monitoring result J{t) by a pure state \^ j) rather 
than by the average state g. This pure states is called 
a conditional state and can be viewed as our state of 
knowledge of the system. However, in systems where in- 
formation about only some of the decay channels can be 
obtained, there is missing information, and it is no longer 
possible to have a pure state description. In this case, we 
can assign a conditional state matrix, qj to represent 
the state of the system under continuous observation of 
the particular decay channels. The evolution equation of 
this conditional state is referred to as a stochastic mas- 
ter equation (SME) [1, 2, 3, 4, 5, 13, 15, 19]. It has the 
property that its average state p(i), 

g{t) = ¥.[ej{t)] (4.1) 

is the solution to the master equation. Here E denotes 
an ensemble average over measurement records J{t) 



In circuit QED, there are four decay channels described 
by the rates Kjn, Kout, 7i and 7^. Only information com- 
ing out of the Kout channel, that is photons transmitted 
through the resonator, are monitored. As a result, we 
must restrict our description to a conditional density ma- 
trix. As will be discussed in the next section, in the 
effective model for the qubit only, this will correspond to 
monitoring the decay channel corresponding to the rate 
rd(t) and not monitoring the 71 and 7^ channels. 

Although, direct detection of the transmitted mi- 
crowave photons is possible [29], here we will consider 
homodyne processing. That is, we will assume that the 
signal coming from the output port of the resonator is 
mixed with a strong local oscillator of phase tuned to 
the signal frequency. Given the homodyne measurement 
result J(t), we can assign to the qubit-resonator system 
the conditional state Qj{t) whose evolution is governed 
by the SME [53] 

Qj{t) =CtotQ.j{t) + i^/Kq[Q4„ Qj]{J{t) - y/Kr}{2I^)t) 
+ ^M[2I^]Qj{t)[J{t) - ^{2I^)t] 

(4.2) 

with £tot given by Eq. (3.1), rj is the measurement effi- 
ciency which we define below, M[c] is the measurement 
superoperator defined as 

M[c]g^ic^{c)t)g/2+g{c-{c)t)/2, (4.3) 

where (c)t = Tr[cgj(t)] and the 0-dependent field compo- 
nents are 21^ = ae^*"^-f ate*"^ and 2Q^ = -me^^'^H-mte*'*. 
The efficiency is rj = Kout^et/'^ with Kout the rate at 
which photons come out of the output port of the res- 
onator and 7/dot is the efficiency at which these pho- 
tons are detected. For the current circuit QED exper- 
iments [7, 8, 9, 24, 29, 30], this can be written as 
'7dct = 1/(A^ + 1) with N the number of noise photons 
added by the amplifier stage. 

For homodyne detection, the measurement record ob- 
served in an experiment can be expressed as 

J(i) = V^(2/^)t+e(t), (4.4) 

where (^ii) is Gaussian white noise and represents the 
photon shot noise. It is formally defined as [40] 

nm\ = (4.5) 

with E denoting an ensemble average over realizations of 
the noise ^(t). 

The measurement term in Eq. (4.2), the one includ- 
ing the superoperator A^[2/0], comprises two parts. We 
will refer to the first one as the homodyne gain and 
the second as the innovation. The homodyne gain is 
^/Kr\M.\2I^gj(i) and collapses the state towards a 
state (eigenstate of I^). The innovation is J — ^J^{2I^)t, 
it pushes the conditional state to a higher or lower 
state depending on whether the current result Jit) is 
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greater or snialler then the average ,jKri{2I^)t- The last 
term in Eq. (4.2), the noisy Hamiltonian term is extra 
non-Heisenberg backaction (it does not come with any 
information gain) caused by the measurement. It is this 
term that stops the conditional state from being driven 
to a state as it supplies random kicks to the con- 
ditional state causing delocalization in J^. 

For heterodyne measurement (where the local oscilla- 
tor (LO) is detuned from the signal frequency), both com- 
ponent and can be measured simultaneously. In 
this case, the SME Eq. (4.2) has extra gain and innova- 
tion terms corresponding to this additional component. 
Heterodyne measurement also lowers the measurement ef- 
ficiency 77dot by a factor of 1/2. Detection of this second 
component does not change the physics in an essential 
way, and we will focus on the simpler case of homodyne 
detection in the remainder of this paper. 

Finally, we note that the SME Eq. (4.2) is known to 
lead to a diffusive- like evolution for the system [3, 13, 41]. 
We will however see in section VI how there can be a 
cross-over from diffusive to jump-like evolution as the 
measurement strength is increased. 

B. Stochastic master equation for the qubit 

In this section, we use the transformation Eq. (3.8) to 
obtain an effective SME for the qubit only. As described 
in Appendix B, this can be done in the limit where 



K 

We note that this limit docs not imply that the informa- 
tion gain about the state of the qubit is small. For the 
system to be in this limit, we simply require that the 
photons come out of the resonator faster than the qubit 
decay rate. Otherwise, the full SME of Eq. (4.2) must be 
used. However, as will be seen from numerical investiga- 
tions, the validity of the model can extend well beyond 
this limit in practice. 

As described in Appendix B, in the limit where 
Eq. (4.7) is vahd, the effective SME for the qubit ob- 
tained using the transformation of Eq. (3.8) and tracing 
over the resonator states is 

pj{t) ^Cpj{t) + yf~(t)XK]p,j(i)(J(t) - yf~(t)(a,)t) 

- *^^^^[^.,pj(t)](J>) - VrWi)(<T.>0- 

(4.8) 

In this expression, rci(t) is the rate at which coherent in- 
formation comes out the resonator and T\^^{t) represents 
extra non-Heisenberg back-action from the measurement. 
These rates are given by 

rci(0 =r/K|/3(0|' cos^icj} - Op) (4.9) 
rba(i) ='nn\(}{t)? sin2(0 " Op), (4.10) 



with the angle Op being defined in Eq. (3.7). 

In Eq. (4.8), J{t) is the processed record coming from 
the resonator and is given by 

j{t) = ^Mz)t+m- (4.11) 

which can be related to the homodyne current by 

J{t) = J{t) + cos(0^ - 4>) (4.12) 

with p{t) = ae(t) + ag{t) and 0^ — arg(/i). 

From the above expressions, by doing homodyne mea- 
surement on the field leaving the resonator wc arc in fact 
doing a weak diffusive measurement of the qubit operator 
Gz with measurement strength rci(t), as well as causing 
extra unitary back-action on the qubit (by the term pro- 
portional to rba(O)- This is clearly seen by rewriting 
Eq. (4.8) in terms of the Bloch vectors 

xj{t) = - {u.,,[t) + v/rba(t)[J(t) - ^V,,{t)zj{t)]}yj{t) 

- {72 + rd(t) + ^T~{i)zj{t)[J - ^T,,{t)zj]}xj{t) 

yj{t) ={^.c{t) + VKJ^[J{t) - Vr~(t)zj{t)]}xj{t) 

- {72 +rd(t) + V'^ciit)zj{t)[J ~ ^T,i{t)zj]}yj{t) 

zj{t) =v/r~(i)[i - zj{tf][j{t) - V^Jijzjit)] 

-7i[zj(t) + l]. 

(4.13) 

From these expression, we see that the effect of the mea- 
surement consists of two parts. First, a nonlinear up- 
date to the state which results in pushing the system 
towards either z ~ ±1 depending on the observed record 
(the terms proportional to rci(t)), and second, an extra 
stochastic component that causes rotations around the 
az axis (the rba(0 terms). We find that the sum of these 
rates is rjV^lt) where 

r,„(i) = K|/3(t)|2 (4.14) 

is the maximum measurement rate. This rate can be 
simply understood by noting that k is the rate at which 
photons leak out the resonator and |/?(t)P is the amount 
of information about the qubit state encoded by these 
photons. That is, this is the maximum amount of infor- 
mation that can be gained from this system via homo- 
dyne monitoring. Looking at Eq. (4.9) this is achieved 
by setting the homodyne phase to 9p and physically cor- 
responds to detecting the quadrature which has the great- 
est separation of the pointer states (see Fig. 2C)). If we 
set (j) to be Op + Tr/2, then rba(i) = V^mit) which re- 
sults in the measurement revealing no information about 
the qubit state and only causing dephasing by qubit fre- 
quency change via photon shot noise. 

When the resonator dynamics have reached steady- 
state, it can be shown that 

rd(i^c») =r„,(t^cx3)/2. (4.15) 

This is illustrated in Fig. 4 where both Ta{t) (solid 
blue line) and r^{t)/2 (dashed red line) are presented 
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FIG. 4: (Color online) r^it) (solid blue line) and rni(f)/2 
(dashed red line) for a measurement drive turned on at the 
time t — 50 ns. The drive frequency is tuned to the bare 
resonator frequency (Ar = 0) and the amplitude is £m/27r — 
\/5. The rest of the parameters are the same as in Fig. 2. 
The inset shows the ratio between the total dephasing and 
the total measurement-induced dephasing. The vertical black 
dashed line is the turn on time for the measurement drive. 
The horizontal dashed line is the steady-state dephasing rate. 



for the situation where a measurement drive tuned to 
the bare resonator frequency is turned on at t = 50 
ns. From this figure, we see that after the resonator 
transients decay away, the dephasing rate is twice 
the maximum measurement rate Fm. Also shown in 
the figure (see inset) is the ratio between the total de- 
phasing of the system, which we define as Fd(s)c?s, 
and the total measurement-induced dephasing, defined 
as /g* Fni(s)ds/2. As it should (when k ^ x), the for- 
mer exceeds, or is at least equal to, the latter. This is 
simply a statement that there is additional information 
about the qubit available in the initial transients of the 
resonator which is not being used when simply monitor- 
ing the homodyne signal. 

A typical trajectory for the conditional state is shown 
in Fig. 5A) [red dashed line is the x component, green 
dashed-doted line is the y component and the solid blue 
line is the z component] for the case of a measurement 
drive turned on at ton = 20 ns and of maximal amplitude 
g.maxy'27r = 5 MHz. The envelope of the drive, shown 
as the green dashed-dotted line in Fig. 5B), is the same 
shape as the one used in Fig. 3. 

Unlike the average evolution, as we monitor the system 
we become more certain of its quantum state. This is rep- 
resented in Fig. 5A) by the z component being stochasti- 
cally pushed towards z — ^l. Note that we are showing 
but one possible trajectory, other realizations will tend to 
localise the z component to -1-1. The fact that our state of 
knowledge about the quantum state is increased by mon- 
itoring is also made apparent by the purity which reaches 
unity [full blue line in panel B)]. Under the average evo- 
lution, see Fig. 3b), we end up with a completely mixed 
state. This is simply due to the fact that the average 
description does not take into account the information 
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FIG. 5; (Color online) A) Numerical solution of the three 
components of the Block vector as obtained from the full ho- 
modyne SME: X red dashed, y green dashed-dotted and z 
solid blue. The measurement frequency and envelope is the 
same as in Fig. 3. We have assumed perfect efficiency rj — 1 
and taken the LO phase to he (f) = 9/3 = ir , corresponding to 
maximal information gain at Ar = 0. B) Conditional purity 
Pj (full blue) and measurement pulse amplitude (dot-dashed 
green). Also shown is one minus the trace distance (dashed 
red), as calculated from the full homodyne SME and the ef- 
fective SME. C) One minus the trace distance minimized over 
the simulation time, 300 ns, for a range of measurement ampli- 
tude and two values of Ti : 7 /is (full red dots) , 100 ns (empty 
blues squares). Here, the trace distance is averaged over 70 
trajectories. 



gain due to the measurement. 

Also shown in this panel B) is the trace distance be- 
tween the conditional state found using the homodyne 
SME Eq. (4.2) and the effective SME Eq. (4.8), in both 
cases using the same parameters as in panel A). The 
agreement between the two results is excellent, which 
shows that when the limit Eq. (4.7) is valid, the effec- 
tive model derived here is indeed a good description. To 
see the break down of the validity of the model, we plot 
in panel C) one minus the trace distance minimized over 
the simulation time, 300 ns, as a function of measurement 
amplitude and for two values of Ti. This is calculated 
and averaged over 70 trajectories. For Ti = 7 /is (full red 
circles), the agreement is very good at all simulated am- 
plitudes. However, for Ti — 100 ns (empty blue squares). 
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the condition Eq. (4.7) is much less valid. As the mea- 
surement amplitude is increased, the approximation that 
we can neglect elements of Eqs. (BIO) and (Bll) in Eq. 
(B9) breaks down and the discrepancy between the effec- 
tive and full models is apparent. In any cases, it is clear 
that the effective model is a very good approximation for 
long Ti and can give a relatively accurate description of 
the system even at low Ti, provided that the measure- 
ment amplitude is not too large. We note that a Ti of 7 
fis was reported in Ref. [8]. The model obtained here is 
therefore useful in realistic experimental settings. 

Moreover, we find that for the small Ti of 100 ns, while 
the discrepancy reported in panel C) for an average over 
many trajectories can be large, when inspecting single 
trajectories we find no notable new features as the mea- 
surement amplitude is increased. The only change can 
be attributed to rcnormalization of the various system 
parameters under the measurement. This suggests that 
the effective model presented here may extend well be- 
yond the limit (4.7), with the only change being that 
the measurement rate must be obtained numerically (or 
experimentally measured) . 



C. Signal-to- noise ratio 

Using the above results, we can define the signal-to- 
noise ratio (SNR) as 



SNR = = yTrnCos^jOp 
71 71 



(4.16) 



As expected, the SNR is maximized by setting the phase 
(j) of the local oscillator to the quadrature containing the 
most information: 4> = Op (see Fig. 2). Since Fm scales 
with the measurement amplitude, and thus with the num- 
ber of photons populating the resonator, this expression 
indicates that arbitrarily large SNRs can be reached in 
principle by simply increasing this amplitude. However, 
as was also pointed out in Ref. [6], at large photon num- 
bers the dispersive Hamiltonian Eq. (2.2) breaks down 
and the results obtained here are simply no longer valid. 

As a result, we must maximize the SNR at a fixed num- 
ber of photons which is assumed to be small enough for 
the dispersive approximation to be valid. The number of 
photons inside the resonator depends on the drive ampli- 
tude Ed, on resonator damping k, resonator-drive detun- 
ing Ar and on the dispersive qubit-resonator coupling x- 
Because of the latter, the actual number of photons can 
depend on the state of the qubit. We thus define Ug and 
rig the number of photons given that the qubit is in the 
ground or excited state [10]. When optimizing, we there- 
fore work at a fixed n = max(ne,ng) which is assumed 
to be well below the dispersive breakdown. Under this 
constraint, and assuming that the resonator has reached 
steady-state, we find that the optimal SNR, for arbitrary 
X and K, occurs at = and is given by 

AnrjKX^ 




FIG. 6: (Color online) The signa-to-noise ratio (SNR) at Ar = 
and for a fixed number of photons in the resonator n = 2. 
The efficiency was taken to be ideal 77 = 1 and the relation 
rate 7i/27r = 1 MHz. The three insets show the phase space 
representation of the resonator state for k = O.Olx, = 
and K = 5x, respectively. The optimal SNR occurs at «: = 2x- 



A plot of this is shown in Fig. 6. From this plot, it is 
clear that the optimal value occurs at k = 2x and is given 
by SNRmax = ■im]xhi- 

Thus, to extract the most information about the state 
of the qubit, we need to tune the measurement drive 
to bare resonator frequency ujr and choose the system 
parameters such that x = k/2. Note that k is chosen at 
fabrication time (and could be tuned in-situ by a change 
of resonator design) and that x = 9^ li'^a — w^) can be 
tuned in-situ by tuning oja- At a drive strength such 
that n ~ ?T,crit [6], the maximum SNR is 77 A/71, which 
for current experiments is approximately 50— 100. Using 
the results of Ref. [23] , this implies that a measurement 
fidelity of at least 95% could be achievable in practice. 
To increase further this value, improvements have to be 
made to Ti, to minimise the amplifier noise and/or tune 
the system to be still further in the dispersive regime. 



V. MEASUREMENT TIME 

In this section, we use the effective SME for the qubit 
to obtain an estimate for the measurement time. That 
is, an estimate for how long it takes us to be confident 
about the state of the qubit. To do so, we will take the 
simplifying assumption that the resonator transient can 
be be ignored (that is Fd = Fni/2). In this situation, the 
simplest definition for a measurement time is to assume 
it is the inverse of the information gain rate: 



— 1/Fci. 

Using this definition, we can write 
imFd < 1/2. 



(5.1) 



(5.2) 



SNR: 



7i(k2/4 + x')' 



(4.17) 



This is the information gain/dephasing 'uncertainty re- 
lation' already discussed by several authors [42, 44, 45? 
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]. The equality is saturated for 77 = 1 and cf) = 9^ and 
corresponds to the quantum hmit discussed in Ref . [44] . 

In the next sections, we define the measurement time 
with respect to three diff'erent measures. As we will see 
however, there is no unique definition of the measurement 
time. 



A. Integrated signal 

For simplicity, let us first discuss the situation where 
71 = and with the qubit initially in either z = ±1. In 
this situation, the measurement time is very well defined 
as the qubit state simply remains in its initial state and 
the task is to estimate that state in the shortest possible 
time. 

We define the measurement signal for an integration 
time t as 

s{t) = ^/f~ I J{t')dt'. (5.3) 
Jo 

This signal has mean and standard deviation 

m = ±tT,, (5.4) 

mt)\ = \l([-<t)-m?) = Vtr:,. (5.5) 

As a result, if we were to measure the system for a time 
t and repeat this measurement several times, we would 
be confident that, to one standard deviation, the value of 
s is itFci ± y/tVci- As a result, to be able to distinguish 
between z = ±1 a reasonable requirement is 

A\s{t)\<tr,i. (5.6) 

The equality defines the measurement time to be = 
l/Fci. In this case, we recover the simple assumption 
discussed above. 

A problem with this approach is that the condition 
which defines distinguishability, and thus the measure- 
ment time, is somewhat arbitrary. Indeed, separation by 
one standard deviation is not the only possible choice 
for distinguishability. Lastly the above discussion only 
answers the question of how much time it takes to distin- 
guish between two predetermine initial states of the qubit 
(under a measurement with Gaussian noise) . A more nat- 
ural question is how long it takes the conditional state to 
collapse to either z = ±1. This can be characterized by 
the conditional variance. 



B. Conditional variance 



It ranges from to 1, with Vj = meaning that we are 
certain that the current value of zj(t) is either ±1 and 
Vj ~ 1 corresponding to complete uncertainty. 

For 7i = 0, Eq. (4.13) can be solved analytically using 
Ito calculus. Doing this yields the solution 

zj{t) = tanh [s{t) + tanh"^ [z(0)]] , (5.8) 

where the signal s(t) is defined in Eq. (5.3). For exam- 
ple, in the situation where z(0) = and arbitrary a;(0) 
and 2/(0) (for example a completely mixed state or an 
x-eigenstate) the conditional variance becomes 

^/(t) =1 -tanh^ [s(i)] . (5.9) 

The conditional variance for a typical trajectory is 
shown in Fig. 7 A) [blue solid line]. Here, we have taken 
Fd/27r = 5 MHz and the qubit to be initially completely 
mixed (x(0) = y{0) ~ z(0) = 0). We see that for this typ- 
ical trajectory, at about 10 ns the measurement is com- 
pleted and the state of the qubit is almost pure. However, 
this is only one typical trajectory and to get a quantita- 
tive idea of the measurement time we must consider the 
ensemble average of the conditional variance. This is the 
red dashed line in Fig. 7 A). 

Taking the ensemble average of Eq. (5.9) gives 

^ _ 1 r°° a:sinh(2a:)exp(~xV2t/) ^^ 

^ " 2\/2^i.3/2 y_oo cosh(2!.) -f- cosh(2x) ^' 

where v = Tdt. Using this, the measurement time tm 
can be defined by requiring that the average conditional 
variance falls below 1/2. Numerically, this is found to 
occur when v = 0.851. Thus, using the conditional vari- 
ance as the measure, the measurement time is found to 
be t„, = 0.851/Fci. 

However, when 71 is non-zero, the conditional variance 
always eventually goes to zero which signifies that we 
are eventually certain that the system is in its ground 
state. This is of course not a useful way to define the 
measurement time. To take into account relaxation, in 
Ref. [23] we studied this question by considering an ef- 
fective classical model for a dispersive quantum non- 
demolition measurement which turns out to be equiva- 
lent to the model presented here when the initial state 
of the qubit is either z(0) = ±1 or a mixture of both 
(i.e., x(0) = y(0) = 0). That is, Eq. (4.13) is the same 
as the Kushner-Stratonovich equation that describes a 
measurement of a two state system with Gaussian noise. 
For completeness and since this is a good approach in the 
situation where relaxation is non-zero, we will therefore 
briefiy review the relevant results of Ref. [23] . 



How fast the ensemble of possible trajectories converge 
to one of the possible values z = ±1 is characterized by 
the conditional variance [46] . This measure is defined as 



VAt) = {al)t - = 1 - z.,{t) 



(5.7) 



C. Initial state fidelity 

In Ref. [23] , to quantify how good a measurement is at 
revealing the initial state of the qubit, we imagine prepar- 
ing the qubit in either the z(0) = ±1 state and then 
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FIG. 7: (Color online) A) Typical trajectories for the con- 
ditional variance (blue solid line) with rd/27r — 5 MHz and 
77 = 1. Average conditional variance (red dashed line) found 
using Eq. (5.10) for the same parameters. Note that the cal- 
culated ensemble average over 1000 typical trajectories falls 
directly over the average conditional variance, and is there- 
fore not explicitly shown. The vertical dash line shows the 
time at which the conditional variance crosses 1/2, tm ~ 13.5 
ns. The green dashed-dotted line is the ensemble average con- 
ditional variance for Ti = 159 ns (7i/27r — 1 MHz) and the 
rest of the parameters are identical to the above. For this 
value of 71, the SNR is 10. B) Average fidelity of the mea- 
surement at estimating an unknown but pure initial state for 
rd/27r = 5 MHz, rj — 1 and 7i/27r — (red dashed line) or 
7i/27r — IMHz (blue dots). The green solid line is a fit using 
Eq. (5.18) with fit parameters Aat = 0.846 (horizontal dashed 
line) and rfit/27r = 14.6 MHz. The two vertical lines corre- 
spond to the times where the fidelity crosses erf(l/\/2). For 
7i/27r = this yields tm ~ 15.9 ns, while for 7i/27r = 1 we 
find tm = 18.5 ns 



generate fictitious records J(t) from that input. Given 
these fictitious measurement records, but now assuming 
ignorance of the initial state, we then ask what was the 
initial state of the qubit. To quantify the efficiency of the 
measurement at revealing the correct input state, we use 
the following fidelity measure F: the number of correct 
assignments minus the number of incorrect assignments 
of the initial state, normalized by the total number of 
assignments. This measure will range from 1, indicating 
that the correct initial state was correctly found for each 
record, to —1, indicating that a wrong assigned was re- 
alized for each record. A fidelity of implies that the 
assignments are completely random. 

The criteria used to make these assignments is the opti- 
mal assignment criteria found in Ref. [23]. This is based 
on an estimate z{t) for the initial state. If this estimate 
is above 0, we assign the initial state as excited and if it 
is below we assign it as ground (for z = 0, we randomly 
chose between excited or ground). This estimate z{t) is 



defined as 

z(t) =P(zo-+l|J,<)-P(zo 



-l|J,t), (5.11) 



where P{zo = ±l\J,t) is the probability at time t that 
the initial state was ±1 given the record J{t). This con- 
ditional probability is given by Bayes theorem 



P{zo\J,t) = 



P{J\zo,t)P{zo) 



(5.12) 



where P{zo) is the initial state probability which we take 
to be 1/2 for both i = — 1 and 1. By introducing a 
fictitious unravelling of 71 (see Refs. [5, 23]), the estimate 
can be rewritten as 



1 - e-2^(*) [1 - a{t)] 
1 + e-2'^« [1 + a{t)] ' 



where 



a(t) =71 / dt^.e-7ife-*)+2^fe) 



(5.13) 



(5.14) 



and s{t) is the integrated signal defined in Eq. (5.3). 

The above integral cannot be solved in general for a 
typical record s{t), but it can be easily determined nu- 
merically. Here, however, to develop a physical under- 
standing of this result, we will consider the 71 = limit 
first. Doing this, we can rewrite z{t) as 



1 



1 + e-2''(*) ' 

and using the above definition of the fidelity 



F(t) = lim — 



E - E + E - E 



z+i>0 z+i<0 z-KO z_i>0 

we obtain the following compact expression 



(5.15) 

1, 
(5.16) 



Fit) 




(5.17) 

This resuh is shown in Fig. 7 B), for rd/27r = 5 MHz, 
as the red dashed line. As can be seen from this figure, 
arbitrarily large fidelities can be achieved by monitoring 
longer. 

We now use this result to define a measurement time. 
A reasonable criterion is to ask for the fidelity to be at 
least erf(l/-\/2) = 0.68 (corresponding to the probability 
of a Gaussian random variable being within one standard 
deviation of the mean). The measurement time is then 
defined as the time when the fidelity reaches this value. 
Using the above results, we find that the measurement 
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time is tm = 1/rci- For 71 = 0, we thus again recover the 
result Eq. (5.1) of the simple model. 

The effect of relaxation on the fidelity is investigated 
numerically. For example, the blue dots in Fig. 7 B) show 
the fidelity obtained numerically for 7i/27r = 1 MHz. 
With all values of Fd tried numerically, we have found 
excellent agreement of the fidelity, including dissipation, 
to the following form 

F(t)=Afitcrf (^y^^ (5.18) 

where Fgt and Afn are fit parameters. For example, the 
green sohd line in Fig. 7 B is for rd/27r = 5 MHz). To 
show how the fidelity F{t) depends on the SNR (which 
in the limit that the resonator transient can be ignored is 
2Fd/7i), numerical simulations and the above described 
fitting procedure were performed for values of rd/27r 
ranging from 1 to 50 MHz. The results of these simu- 
lations are shown in Fig. 8. In particular, the inset of 
this figure shows the long time fidelity, A^t- From this 
result, we see that the fidehty is below erf(l/v^) until a 
SNR of 3.16 is reached (this corresponds to the vertical 
dashed line). Below this value, using the above measure- 
ment time criteria, we must conclude that the measure- 
ment cannot be completed. However, for all SNR above 
3.16 it is possible to find a time where the fidelity crosses 
erf(l/-\/2)- This time is shown in Fig. 8 as the blue solid 
line. We see that this time is always longer then the ideal 
limit (red dashed line). However, as the SNR increases, 
the measurement time approaches the result obtained in 
the 71 = case (l/Fci). 

VI. EMERGENCE OF QUANTUM JUMPS 

By increasing the SNR (e.g. increasing the measure- 
ment drive amplitude), the conditional state will tend to 
be localized to cither the excited or ground state. This is 
due to the term vT^(l — Zj) in Eq. (4.13). However, in 
the presence of relaxation 71 7^ 0, on average the condi- 
tional state must be localized on the ground state in the 
long time limit. As a result, when initializing the system 
in the excited state, there are two competing processes. 
One is trying to lower the conditional state to 2 = — 1 
(71) and the second is trying to hold it in at z = -1-1 (Fci). 
The net result is that as the measurement strength be- 
comes stronger, the SME will transform from a diffusive 
stochastic process which stochastically lowers itself to the 
ground state to one which jumps at a well defined time. 
That is, as the measurement becomes stronger, quantum 
jumps due to the decay channel 71 will be revealed in the 
signal J{t). 

The crossover of the conditional state from a diffu- 
sive to jump-like SME is numerically shown in Fig. 9 
for measurement drive strengths ed/27r = {0, 0.2, 0.9, 20} 
MHz, respectively. From this figure, we see that when 
ed/27r = (dashed red in panel A), as expected, the 
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FIG. 8: (Color online) Measurement time (full blue) as 
extracted from the requirement that the fidelity reaches 
erf(l/^) as a function of the SNR (2rd/7i). This is plotted 
for 7i/27r = 1 MHz. The vertical dashed line corresponds to 
SNR — 3.16. This defines the minimal SNR for which we can 
define a measurement time. This is also presented in the inset 
which shows the long time fidelity as a function SNR, with 
the vertical dashed line being the SNR where the long time 
fidelity crosses erf(l/\/2). The red dashed line in the main 
panel corresponds to tm ~ l/Fci the measurement time for 
71 = 0. At low SNR, it always underestimate the measure- 
ment time. 



state simply decays to the ground state at a rate 71 as 
this is equivalent to gaining no information about the 
qubit state. For ed/^n ~ 0.2 MHz (full blue in panel 
A) it stochastically follows this decay curve. However, at 
ed/27r = 0.9 MHz (full blue in panel B), there appears 
a sharp transition from the excited state to the ground 
state. As the measurement amplitude is further increased 
(ed/27r = 20 MHz, red dashed in panel B), this transition 
becomes sharper and happens at a well defined time. 

Experimentally, one does not have direct access to zj (t) 
but reconstructs this quantity from the measured current 
J(t). As a result, a more striking demonstration of quan- 
tum jumps would be to see these jumps directly in J{t). 
However, due to white noise fluctuations ^(t), without 
first averaging the signal over some finite integration time 
the jumps are not visible in most trajectories. As a result, 
we consider the quantity 




where dt is an adjustable integration parameter. At large 
St, the fluctuations are well suppressed by the averaging, 
but at the expense of smoothing out the jump. For exam- 
ple. Fig. 10 shows the integrated current corresponding 
to the two measurement amplitudes used in Fig. 9 B). 
Panel A) of Fig. 10 corresponds to ed/^n = 0.9 MHz and 
St = 48 ns, while panel B) corresponds to ed/27r = 20 
MHz and St = 32 ns. From this flgure, we see that for 
Ed/^TT = 0.9 MHz one cannot resolve the jump from the 
integrated signal. However, by increasing the drive to 
ed/27r = 20 MHz the jump becomes clearly visible. 
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FIG. 9: (Color online) Time evolution of the z component of 
the conditional state for a qubit initially in the excited state 
and x/27r = 5 MHz, k/2tt = 10 MHz, = 0, r? = 1 and 
Ti — 7 fis. A) Measurement drive strength of ed/2n = MHz 
(red dashed) and 0.2 MHz (solid blue). B) 6^/2% = 0.9 MHz 
(solid blue) and 20 MHz (dashed red). 



FIG. 10: (Color online) Integrated signal S{t) as a function 
of time. A) Measurement drive strength Ed/^n = 0.9 MHz 
and integration steps 5t = 48 ns B) ed/2n = 20 MHz and 
St = 32 ns. The vertical dashed line indicates the time where 
the conditional state crosses in Fig. 9. The two horizontal 
dashed lines are the expected results corresponding to ^ = ±1 
constant. 



Thus wc sec that if wc wish to sec the jump in the 
integrated signal we need a much higher signal to noise 
ratio. This is because by definition the conditional state 
is the best way we can process the record. Wc note that 
it should be possible to optimize the processing of the 
current to distinguish the jump at lower SNR by using a 
more complicated filter then time averaging. For exam- 
ple, one could easy envisage that there exists an optimal 
integration time St — topt, or kernel k{t), such that the 
integrated signal, S{t) = J^^k{t,s)J{s)dt, would do al- 
most as good as the conditional state (see Ref. [23] where 
a similar question was investigated). 

Finally, an interesting situation occurs if the 71 channel 
can be directed (corresponding to the fluorescence of the 
"atom"). When we register that the qubit has jumped, 
we know that the energy has been dissipated into the 
71 channel. If this fluorescence is directed to an addi- 
tional low-Q resonator, this can be used as a single mi- 
crowave photon source with a well-defined emission time. 
Inversely, if the qubit is prepared in the ground state this 
second resonator can be used as a single photon detector. 
When a jumped of the qubit is registered, we know that 
a photon was present in the second cavity. The efficiency 
of this detector will be limited by the Zeno effect [54] . 



VII. WEAK DRIVING OF THE QUBIT - ZENO 
EFFECT 

In the previous sections, we considered the dynamics 
of the qubit under measurement. In this section, we will 



include the possibility of having a control drive. As ex- 
plained in Sec. II, this is done in practice by adding a 
drive on the input port of the resonator at the qubit 
transition frequency. Our goal is to use the transforma- 
tion Eq. (3.8) to again obtain an effective SME for the 
qubit only. However, as shown in appendix C, in the 
transformed frame, the drive on the qubit induces tran- 
sitions between the different (transformed) states of the 
resonator. As a result, it is not possible to integrate out 
exactly the resonator from the effective description of the 
qubit. 

However, in the small P limit, we find we find that 
qubit evolution is well described by the SME 

Pjit) W=c,Pj{t)] + [70 + Td{t)] V[a,]pj{t)/2 

^r^[a^,pj{t)]{J{t)-^,{a.)^) 

+ ^/f^,M[G,]pj{t){J{t) - Vr^(a,)t) 
+ 7il?[a_]pj(0 

(7.1) 

where VLn is the Rabi frequency of the control drive tuned 
to the qubit frequency, Wac- The competition between 
Rabi flopping induced by the control drive and the extra 
dephasing due to the measurement drive should enable us 
to observe the Zeno effect in this system. Indeed, as the 
measurement rate Fci is increased, the qubit dynamics 
will switch from Rabi oscillations to jump-like behavior. 
In the limit of very large measurement rates, all dynamics 
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will disappear and the qubit will remain fixed in the zj = 
±1 state (if 7i is zero). 

To see this more clearly, the above SME can be written 
in terms of the components of the Bloch vectors as 

ijW = - {72 + Tdit) + Vf~,zj{t)[J{t) - x/f~,zj{t)]}xj{t) 
yjit) = - QRZjit) ~ {j2 +Td{t) 

+ V^iZj{t)[Jit) - Vr~zj(i)]}yj(t) 
zjit) =nRyj{t) - -fi[zjit) + 1] 

+ v^[l - zj{tf][Jit) - Vf~zj(i)]. 

(7.2) 

Here we have assumed that (p — 9p such that Fba = 0. In 
the limit that the measurement induced dephasing rate 
Fd is much larger than all other rates, we can set xj{t) = 
yj{t) = to arrive at 

2j(i) = - 27jump2,7(t) - 7i(z,/(t) + 1) 

+ ^,{l-zj{tf){J{t)-^,zj{t)). 

As expected, in this expression, the jump rate is given 
by 7jump = ^^fl/2(72 +rd)- From this expression, we can 
indeed expect Zeno-type dynamics for the qubit. Note 
that the stochastic term \/T^iZj{t)[J{t) — \/TdiZj{t)] has 
been ignored in the denominator of the jump rate. 

As discussed in Sec. VI, as the measurement becomes 
stronger (i.e. Fci increases) the conditional state will lock 
itself into either the excited or ground state {z = ±1). 
Because of the control drive, whose amplitude enters in 
the jump rate 7jump, random fluctuations will result in 
the system jumping between these two localized states. 
However, unlike the situation presented in the previous 
section, the rate of the these jumps also depends on the 
measurement strength through Fci. As a result, in the 
limit that the measurement SNR is large, the jumps will 
disappear and the qubit will remain fixed in either z = 
±1. This is the Zeno effect [25, 26, 27, 28]. 

To illustrate this, we have numerically integrated the 
full homodyne SME (4.2), with the addition of a qubit 
drive, for 2 different measurement drives ed/^ir = 0.9 
MHz and 20 MHz, and a Rabi frequency of flii/2TT = 
2.5 MHz. The results are shown in Fig. 11. From this 
figure, we see that for the low measurement amplitude 
(pannel A), measurement only causes small amplitude 
noise on the Rabi oscillations. However, for the larger 
amplitude (panel B) the conditional dynamics is jump- 
like. We note that the rate of the jumps between the 
ground and excited state is, as expected from the above 
discussion, slower then the rate Qr at which the Rabi 
drive would coherently drive the system. Note that the 
X-axis scale is not the same in pannel A and B(C). 

These jumps can be observed directly in the signal 
J{t). This is illustrated in Fig. 11 C) for the case of 
the larger measurement amplitude. As in the previous 
section, we are showing here the integrated signal, de- 
fined in Eq. (6.1), and with a integration time St ~ 32 
ns. The integrated signal clearly reveals the jumps. 
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FIG. 11: (Color online) Time evolution of the conditional 
state for the full homodyne SME for a qubit initially in the 
ground state: x red dashed, y green dashed-dotted and z solid 
blue. The system parameters are x/27r = 5 MHz, k/2tt = 10 
MHz, = 0, r? = 1, 71 = 70 = and Qn/^TT = 2.5 MHz. 
A) Measurement drive amplitude of ed/2n = 0.9 MHz. B) 
ed/27r — 20 MHz. C) Integrated signal with St = 32 ns and 
for a measurement amplitude of £d/27r = 20 MHz. 



We note that the numerical results presented here were 
obtained using the full homodyne SME, and not the sim- 
plified SME (7.1). This is because while this simplified 
equation contains the essential features, it is not strictly 
valid in the limit of large mean photon number. How- 
ever, by inspecting individual trajectories, we find that 
trajectories obtained from the full model display the same 
qualitative behavior as those obtained from the effective 
model. The difference between the two types of trajec- 
tories appear to be only small changes in renormalized 
system parameters. 

To explore the difference between these two models in 
more detail, we simulated both models as a function of 
the measurement drive strength Sd- The average jump 
time extracted for both cases is shown in Fig. 12. The 
simple model (green filled diamonds) agrees with the pre- 
dicted jump time l/7jump for large £d (blue solid line). In 
the low measurement drive limit, it agrees with the time 
set by the Rabi drive (red dashed line). However, the full 
homodyne SME calculations (red open diamonds) do not 
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FIG. 12; (Color online) Average jump time as a function of 
the measurement drive strength Ed- The jump time is defined 
as the time for the system state to go from z — ±0.95 to 
z — =f0.95. The red dashed line corresponds to the Rabi 
period tt/JIr of 200 ns for a Rabi drive of flii/2n = 2.5 MHz. 
The blue solid line corresponds to the expected mean jump 
time expected from the Zeno effect, 7jump = f2|;/2(72 + Fd). 
This agrees very well with the numerically simulations of Eq. 
(7.1) (green filled diamonds) where the mean jump time was 
extracted from a total simulation time of 16 ms. The red 
open diamonds are the mean jump time as extracted from 
the full numerical homodyne SME for a total simulation time 
of about 15 ms. 



agree with this expectation in the large drive case. Al- 
though the mean jump time does increase with the mea- 
surement drive as expected by the Zeno effect, its effect 
is not as large as that predicted by the simple model. . 



VIII. CONCLUSION 

We have investigated the measurement of a super- 
conducting qubit using a dispersively coupled resonator. 
From the fully quantum mechanical model for the res- 
onator and qubit dispersive interaction, we have found an 
effective master equation for the qubit (only) under the 
presence of a measurement drive on the resonator. This 
was done by moving to a frame which takes into account 
the entanglement between the resonator and the qubit. 
With respect to the bare qubit-resonator master equa- 
tion, the effective qubit master equation contains an addi- 
tional, time-dependent decay channel which has the form 
of a dcphasing process and which depends on the num- 
ber of photons populating the resonator. This is referred 
to as measurement-induced dcphasing [10, 42, 44, 45? 
]. In addition to the extra dcphasing channel, the pho- 
ton population of the resonator is also responsible for an 
ac-Stark shift of the qubit transition frequency. As was 
shown in Refs. [10, 38], and later reported experimentally 
in Ref. [24] , these effects lead to number splitting of the 
qubit spectrum in the limit where x ^ We thus have 
obtained a very simple to solve (2 dimensional Hilbert 
space) effective model for the qubit dynamics that con- 



tains the essential physics. We believe that this is an 
interesting analytical, and numerical, tool to further ex- 
plore the dispersive regime. 

We then focused on the situation where one is con- 
tinuously monitoring the voltage at the output port of 
the resonator. In that situation, the conditional state 
of the combined system evolves according to the homo- 
dyne SME [13]. Again moving to a frame that takes 
into account the qubit-resonator entanglement, we then 
obtained an effective SME for the qubit that is valid 
in the limit 71 ^ k. The SME corresponds to a weak 
measurement of the qubit operator cr^. As a result, the 
measurement stochastically projects the conditional state 
into eigenstates of cr^. This is similar to the result that 
was obtained in Refs. [14, 15, 16, 17, 18, 19, 20, 21] for 
a double quantum dot qubit monitored continuously by 
a quantum point contact (although the results are pre- 
sented in very different forms in all of these references). 

From these results, we then showed how, in the limit 
that the measurement becomes strong, there is a cross- 
over in the qualitative behavior of the trajectories from 
diffusive to jump-like. This process relies on the Ti de- 
cay channel for the qubit to jump from the excited to 
the ground state. As a result, we expect to see quantum 
jumps, with Poisson statistics, in the measured homo- 
dyne current. 

Moreover, in presence of a Rabi drive on the qubit, 
we again saw that the trajectories evolve from diffusive 
to jump-like. In this situation and when the measure- 
ment drive is not too strong, the jump rate is given 
ri|./2(72 -I- Fd). This is due to a competition between 
the measurement drive that tends to lock the conditional 
state on one of the eigenstates of cr^ and the Rabi drive 
that tends to flip the qubit. This is nothing but the 
Zeno effect. However, as the measurement gets stronger, 
the simple dispersive Hamiltonian approximation breaks 
down and the effective model fails to predict the correct 
mean jump time for the qubit. We found that the form of 
the trajectories are still consistant with the model, only 
the parameters appear to be renormalized by the stronger 
measurement drive. An interpretation of this discrep- 
ancy is that the resonator response is not instantaneous. 
When the qubit jumps, the resonator state will lag be- 
hind and as a result the effective separation in phase 
space between the two coherent states of the resonator 
corresponding to the two qubit states will be reduced. In 
turn, this can lead to a reduction of the measurement 
induced-dephasing. 

Overall, the model obtained here can be used to ob- 
tain qualitative information about the qubit dynamics 
in the presence of dissipation, measurement and control. 
We also gave a numerical prescription to go beyond this 
simple model. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE QUBIT MASTER EQUATION 

In this appendix we show how to obtain the effective 
qubit master equation Eq. (3.10) from the full dispersive 
qubit-resonator master equation Eq. (3.1). For this pur- 
pose, we introduce the the transformation 



P = IleD{ae) + ngD{a 



9/' 



(Al) 



where Ile,g are the projectors on the qubit excited and 
ground state and D{a) is the field displacement opera- 
tor [47]. 

In this frame, the state of the combined qubit- 
resonator system = gF can be written in the energy 
basis as 

oo 

n,m— i,j—e,g 

Using this expression for the system's density matrix in 
the transformed frame, it is simple to express the qubit's 
reduced density matrix in laboratory frame as 

p = Trres.[Pe''P^] 



(A3) 



Y^[\n,m,m,n\e){g\ + A;,_„_„_,„ jg) (e|] 



with 



A 



n,m,p,q — Qn,m,e,g 



dp,q exp [-2lm[aga*]] (A4) 



and dp q is the matrix element of the displacement oper- 
ator in the photon number basis 



dp,, = {p\D{/3)\q). 



(A5) 



In the remainder of this appendix, we show how by 
using the transformed master equation for it is pos- 
sible to obtain the coefiicients entering in the expres- 
sion Eq. (A3) for p. 



In the displaced frame, the state matrix g^ obeys the 
master equation 



P^Pg^ - g^P^P 



(A6) 



with ~ P^^OP the transformed operators and where 
the effective Hamiltonian iJoff is given by Eq. (2.2). Using 
the standard results 



D'^{a)aD{a) = a + a 
D\a)a)D{a) ^ a) + a* 



(A7) 



we obtain 



a = a + IIq 
{a^af = a^a + a^Jla + all* + jagl^IIg -f- jael^IIe 

where we have defined 

= agXlg + ajle- (A9) 

For the qubit operators, we obtain 



p 

cr, = a. 



a-D\ag)D{ae). 



(AlO) 



Using these results, we then have for the transformed 
Hamiltonian 

^^cff = + ^ [ed{a^ + n;) + h.c] 

-I- h{/\r + x(^z) [a^a + fltn^ + an; + |n„|2) 

(All) 

We now turn to the damping superoperators. For the 
field damping; we find 



1, 



V[a'-]g'- ^V[a]g'' - -[a^n„ - aH^, g^] 

P* r P 1 , r Pi t 

[az, g \, 



-V[a,]g'^ -i- 



where 



and 



Fm = k\. 



P = Ue- ag 



(A12) 
(A13) 
(A14) 



In this frame, damping of the field induces extra dephas- 
ing of the qubit at a rate T^. Note this is not the mea- 
surement induced dephasing rate given in Ref. [10], as the 
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present expressions arc in the transformed frame, not the 
laboratory frame. Transforming the qubit dissipation su- 
peroperators, we obtain 



transition frequency as 



(A15) 



V[a^]Q^ = V[a-D\ag)D{a,)]Q^ . (A16) 

The last terms to take into account are those in the 
last line of Eq. (A6). For these terms, we use the fact 
that 

D{a) = [{aa^ - a* a) + [a* a - aa*)/2\ D{a) (A17) 
to obtain 

iino} - n* a, 



(A18) 



Choosing dg and dg in the above expression as in 
Eq. (3.5); we finally obtain the transformed master equa- 
tion 



B^{l3e*, + P*ed)l2 + Klu,[agal]. 



(A20) 



In this transformed frame, the last two terms of Eq. (A19) 
appear because the resonator acts as a non-Markovian 
bath for the qubit [48, 49]. Moreover, we note that the 
term proportional to 71 in Eq. (A19) can be expressed as 



V[a-D\ag)D{a,)]Q^ 

= a^D{p)Q^a+D\l3) - a+a^g^ /2 - g^a+a./2. 

(A21) 

Upon taking the trace over the resonator space, this ex- 
pression will simply take the form of regular Markovian 
71 damping on the qubit. This is one of the features of 
the transformed master equation that will allow us below 
to obtain a simple effective master equation for the qubit 
only. 



.p .<^a + -B r p 

Q 7, k^'fi' 



i [(Ar + X(Jz)a^a, g^] 
-iiV[a^D\ag)D{a,)]g^ + + Tj2]V[a,]g^ /2 
KD[a\g^ +'^[az,g^]u, -r — u. 

(A19) 



where we have defined the time-dependent shifted qubit 



2. Moving back to the laboratory frame 

Using Eq. (A2) and Eq. (A19), we get the following set 
of coupled differential equations for the coefficients of the 
laboratory frame reduced qubit density matrix Eq. (A3) 



J 



■ p 

Qn. 



m.e.e 



-iK{n - rn) - ix{n - m) - 71 - K{n + m)/2]g^,„ + K£'^+i,,„+i,e,e V (" + !)("* + 1) (A22) 

Sn,m.g.g =[-iAr(n - m) + ix{n - m) - K{n + m)/2]e^,„_g g + 7i XI ep,g,e,edn,pd*„,^g 

p,q (A23) 



^n,m,p,q =(>n,m,e,gdp,qG ilmlagoj _ ii9j(Im[Q:ga*])A„_m,p,i} + l3y^Xn,m,p-l,q — P* y^Xn,m,p,q-l — -^9tiP P* ) '^n,m,p,q 



(A24) 

~i{u}a + B) - iAr{n - m) - ix{n + m) - 71/2 - 70 - Ed - K{n + m) /2]\n^m,p,q 



+ KXn+l,m+l,p,q\/ (n + l){m + 1) + KfWm + lXn,m + l.p,q - Kf3* VnTlXn + l,rn,p,q + P y/pX„^rn,p-l,q (A25) 
P ^/QXn^m.p.q—l 

I 

where (A24) was obtained by differentiating (A4), with To reconstruct the effective master equation from these 
respect to time and amplitudes, we note that the time-derivative of the |e)(e| 

component of Eq. (A3) can be written as 

Ed = 2xImKa:] (A26) 

and 

B = 2xRc[aga^]. (A27) V V 
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For the ground state component we have 



SME Eq. (4.2). This is 



(A29) 



where in the last equahty we have used the identity 



n,q — 'Jpq- 



(A30) 



n=0 



The off diagonal terms are the hardest to deal with, but 
because in the transformed frame, the photon population 
is initially zero and there is no mechanism to populate 
^n,m.p.q tcims with {n,TO,p, (/} > we have 

Pe,g = Ao,0,0,0 = [-i(l^a + B) - Ji/2 - 70 - Td]pe.g- 

(A31) 

The above expressions correspond to the following ef- 
fective master equation for the qubit 



where Wac = + B. 



(A32) 



gj =Ctotgj + y/Kr]M[2I^]gjJ + i^^[Q^, gj]J. (Bl) 



where the bar here is used to signify that the state is not 
normalized and the linear measurement superoperator is 



M[c]g ^ cg/2 + gc/2. 



Moving to the frame defined by Eq. (3.8) gives 



g^ =Aot£'j + ^/Krj \(3\ cos(6'^ - (t>)M[a^]g^ 



(B2) 



\p\cos[ 



)gj + ae ^ gj + gja'e 



J (B3) 



K-q\l3\ sin(6'/3 - (j)) _p 

— 7^ Wz^g \J- 



APPENDIX B: DERIVATION OF THE 
EFFECTIVE QUBIT STOCHASTIC MASTER 
EQUATION 

To derive the effective qubit stochastic master equation 
we start by using the linear form [5, 50, 51] of the full 



In the same way as Eqs. (A22) - (A25) were obtained, we 
can write equations of motion for the coefficients gn^m,e,ei 
gn.m,g,g \i.m,p,q of thc decompositiou of the condi- 
tional density matrix in the energy basis. The first 
term of Eq. (B3) yields the same terms as in Eqs. (A22) 
- (A25). As a result, we get 



en,m,e,e =(A22) + ^[Vn + le "^gn+l,m,e,e + Vm + le*"^£»^^,„+i_g_g + C0s(6'^ - (l))gn,m,e,e 
+ \P\cOa{dp - (t))gn.„,e.eU 



cos 



-n,m,g,g 



£'n,™,3,3 =(A23) + yH^Wn+le "^gn+l,m^g^g + Vm+^e"^gn,m+l,g,g + \^^\ 
- \P\cOa{dp - 0)Cm,3,g]-^ 

Xn,ra,p,q =(A25) + V^iV^er^'^' g^_,,^^^,^gdp,qer^''-^^'^^<^ + V^^e^*^^^„+i,,^gdp,,e-'i-["«<l 

+ \fl\ C0s{9fj, - (j))Xn,m,p,q + sin{9p ~ 4>)Xn,m,p,q]J- 



(B4) 
(B5) 
(B6) 



In thc above expressions, the equation numbers refer to 
thc RHS of the corresponding expressions. They are the 



contribution of the Lindblad term C 



totgj ■ 



By noting that these matrix elements are not coupled 
in the qubit basis, we can use the same arguments that 
were used to obtain Eq. (A31). Doing this, we find for 



thc stochastic parts 



£'0,0,e,e 



/K?7 



\p\cos{ 



Pe.g — '^O, 0,0,0 



I cos{9p - (j>) 



KT] cos{d^ - 0) + i\f3\ sin(6'/3 - (f>) 



X Pe,eJ 

(B7) 

X Pe,gJ- 

(B8) 
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For the ground state amplitude we have 

Pa, a ~ ^ y Sn,n,g,g 
n 

=yK^p/^| cos(6'^ - 0) - \P\ cos{ep - (f)) 
where we have defined 



(B9) 



-p 



n. 



(BIO) 
(Bll) 



The amphtude gf^ obeys the equations of motion 



+ cos(6'^ - 0) - |/3| cos(6'/3 - (^) 



E 



(B12) 



and we find a similar equation for dtgfr- From this result, 
we find that provided that 



e = < 1, 

K 



(B13) 



the efi^ect of gf i and p^'r on g "^iH be very small. In this 
situation, it is possible to construct a laboratory frame 
linear SME for the qubit only: 



pj =Cpj + ^/Krj\(3\ cos(6'a - (l))M[az]pjJ 

,y^\f3\sm{e0 - 0) _ 
+ » 2 ['^^^PJ\J 

+ y/Kr]\p\ cos{e^ - 0)pj J- 



(B14) 



Using Eq. (4.12) and normalizing the above SME gives 
Eq. (4.8) with measurement record statistics given by 
Eq. (4.11). 



APPENDIX C: INCLUDING COHERENT 
CONTROL OF THE QUBIT 

A drive on the input port of the resonator at a fre- 
quency ujc close to the qubit transition frequency can be 
used to coherently control the state of the qubit. As 
shown in Refs. [6, 31], this can be represented by the 
term 



(CI) 



in the qubit-resonator Hamiltonian Eq. (2.2). In this 
expression, flu is the Rabi amplitude which depends on 
the drive amplitude and its detuning to the resonator 
frequency. 

To see how this changes the effective qubit master equa- 
tions obtained in the previous appendices, we apply the 
transformation of Eq. (3.8). In a frame rotating at the 
frequency frequency ujc, we find 

In this transformed frame, it is apparent that when the 
qubit flips, it must 'drag' the photon field populating the 
resonator. In this situation, the equation of motion for 
the coefficients of Eq. (A3) arc modified in the following 
way 



• p 



— (A22) + i—— ^ ^ [An,p,p,Tn ~ ^ 



(C3) 



(A23) - ^[Xp,n,m,p - >^p,n,m,p\ (C4) 



iA25) -^^y\g^ 



l,m,g,g'^l,' 
I 



(C5) 



where the equation numbers refer to the RHS of the cor- 
responding expressions. These additional terms mix the 
various coefficients of the decomposition of p and obtain- 
ing an effective master equation for the qubit only is no 
longer possible exactly. This can nevertheless be done in 
the small measurement amplitude limit. 

The case where a control and a measurement drive are 
acting on the system simultaneously is most interesting 
in the case where the measurement drive is of small am- 
plitude. For example, in Ref. [8], a weak continuous mea- 
surement drive corresponding to about one photon was 
used to monitor Rabi oscillations. When the measure- 
ment drive amplitude is increased, the qubit is dcphascd 
faster and coherent control is realized with less fidelity. 

In the weak measurement limit (3^0, the matrix 
clement of the displacement operator dn_rn in the photon 
number basis can be approximated to 



/min(p,gf) 



X L 




min(p,g) vl/^l / 



-'^'^/%/^^^[l-(max(p,,) + l/2)|/5n 
y max(p, q)\ 



f3P-\ p>q 
q>p 



(C6) 



where L'^{x) is an associated Laguerre polynomial. To 
lowest order in /?, this reduces dp ^ ~ Sp,q- Making this 
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replacement in the above expressions for the amplitudes, 
we find in this limit the following effective qubit master 
equation 

P = -l ^ Wz,p\''l—[(Tx,p\ ^^^^ 

+ iiV[G^]p + [7^ + Td] V[a,]pl2. 



Unsurprisingly, in the small /3 limit, the only effect of a 
control drive is to flip the qubit at the Rabi frequency ^r. 
If we were to include the next order m. dp^q we would not 
be able to obtain an equation for just the qubit, effects 
such as sidebands would be observed. 
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